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A detailed Linear Combination of Atomic Orbitals (LCAO) tight-binding model is de- 
veloped for the layered High- Temperature Superconductor (HTSC) cuprates. The band 
structure of these materials is described using a u-band Hamiltonian employing Cu4s, 
Cu3d a .2_j ; 2 , 02p x and 02p y atomic orbitals. The Fermi level and the shape of the 
resulting Fermi surface are fitted to recent Angle Resolved Photon Emission Spec- 
troscopy (ARPES) data to realistically determine the dispersion in the conduction band. 
Electron-electron interactions and, ultimately, Cooper pairing is obtained by introduc- 
ing a Heitler-London, two-electron exchange between adjacent orbitals within the Cu02 
plane. Finally, using the LCAO wavefunctions determined by the band structure fit, the 
Bardeen-Cooper-Schrieffer (BCS) type kernel is derived for interatomic exchange. 
Keywords: layered superconductors 

1. Introduction 

Since the first discovery in the high-temperature superconducting cuprates 

have attracted great attention.Em There is a wealth of experimental data about this 
unusual family of materials which we shall not attempt to review here. Rather 
we will mention only two lines of inquiry: the Angle Resolved Photoemission 
Spectroscopy (ARPES) or Angle Resolved Ultraviolet Photoclectron Spectroscopy 
(ARUPS) experiments to determine the normal state band structure, and the 
ARPES measurements of the gap in the superconducting state. From these ex- 
periments there are quantitative results for the shape of the Fermi surface and the 
angular dependence of the superconducting gap around this surface, both of which 
can be compared in detail to any proposed theoretical model. 

The quantity of theoretical work on the cuprates is similarly large, but can be 
roughly divided into two camps: one has made considerable effort to realistically 
model the band structure of the materials, from first principles where possible, 
whilst the other has concentrated on correlation effects by using idealized "toy" 
models of the Hubbard and t-J types. Whilst both of these camps have made 
considerable progress in their own fields, there has been surprisingly little overlap 
between the two. In particular, the band structures employed in the studies of 
correlations are rarely more sophisticated than the inclusion of nearest- or next- 
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nearest-neighbour hoppings on a square lattice, and bear little resemblance to either 
the measured or predicted band structures. 

It is the purpose of this paper to analyze microscopic interatomic two-electron 
exchange processes which could be useful to understand the nature and origin of the 
superconducting correlations in the cuprate materials. As it is now almost univer- 
sally accepted that it is the physics of the CuC"2 planes, the common feature of all 
the cuprates, which gives the materials such extraordinary properties our efforts will 
focus on the properties of these planes and we will omit the effects of extraneous, 
material-dependent details such as the presence of CuO chains, orthorhombic dis- 
tortions, double planes, etc. Having devised a model we will fit recent experimental 
data to the single-particle sector, so that it accurately describes the band structure, 
and then we will proceed to derive an analytic expression for the BCS-type ker- 
nel, or pairing potential, V(p,p'). We will leave the solution of the resulting gap 
equation, and comparison to experiment, to a later date. 

2. Model 

In this work we attempt to understand the physics of the CuC-2 plane using a 
relatively simple, idealized, tight-binding model, within the Linear Combination of 
Atomic Orbitals (LCAO) approximation. To this end we consider only the four 
atomic orbitals coming from the Cu3d x 2_ y 2 , Cu4s and 02p CT states. Denoting the 
positions of the Cu, O x and Oj, atoms by Rc u , R-z and H y respectively, the in- 
plane lattice constant by ao and the unit cell index by n = (n x ,n y ), the LCAO 
wavefunction reads: 

V'LCAoO) = ^2 [ D nai>Gu3d x 2_ y 2 - nfl - Rc«) + S na 1pCui s (r - na - R Cu ) (1) 

+X na tp 2 P!C - na - R x ) + Y na ip 02py (r - na - R y ) , 

where $„„ = {D ncn S na , X nai Y na ) is the vector of amplitudes for the nth unit cell. 
With each of these quantities is associated a Fermion operator satisfying the usual 
anticommutation relations {X na , X^} = 5 nm S a p, {A nQ , Y^} = 0, etc. We use 
these operators to analyze exchange processes in the Cu02 planes. 

This model can be conveniently split into two parts; a non-interacting part 
describing the single-electron band structure and an interacting part describing the 
effects of electron-electron interactions. We consider each of these in turn below. 

2.1. Tight-binding band structure 

When the differential overlap of the orbitals is neglected, the non-interacting, single- 
electron Hamiltonian can be parameterized by three atomic energies, e s , and 
e p , and three hopping integrals, t sp , t p d and t pp , corresponding to the transitions 
Cu4s <-» 02p a , 02p a <-> Cu3d x 2_ y 2 and 02^ <-> 02^ respectively. In second- 
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quantization notation, the resulting Bloch-Hiickel-like Hamiltoniani reads 

Hb,U = ^ S -D„a Cd£) nQ — t pd ( — X na + X( n;E _ lny)a — Y na + Y(n x ,n v -l)a) (2) 



* not 



j5*na ^-sp ( ^ nrv H~ -^(n x — 1,72.^)0: ^ nrv H~ ^(rij. ,n y — l)a 

£pp(^na ^~ ^(nj + ljiyja "~ [n Xi n y — l)a ^(n^+l,^ — l)a) 
Cp^na tpd[ D na + ^(n x ,n y + l)a) ^sp ( $na 4~ $(n x ,n y -\~l)a 

^pp(-^na ^(n^ — l,n y )a ,n y + l)a ~t~ -^(n x -l.n y + l)a 



This Hamiltonian is diagonalized by introducing the Bloch states $ pa = (D pa , S pa , X pa , Y pa ). 
The elements of these states are defined by the relation 



( ^ na \ 

Sua 

V Y na J 



= 1 Wp- 



•'pa 



(3) 



where TV is the number of unit cells, and the two phases are <p a = h(Px — t) and 
0b = 2 (Py ~ 1") • The Hamiltonian can now be expressed in the more compact form 



H B H = Y,* Pa H p V pa , 



(4) 



where the Hamiltonian matrix is 



( e d t p dS x 

tpdSx v S pS x £p 
V ^pd^y t'spSy ^pp^x^y 



tpd^y \ 



isp$y 
~tppS x Sy 



(5) 



Here, for the^ake of convenience, we have adopted the shorthand notation of An- 
dersen et al.JaB s x — 2sin(p x /2), s y = 2sin(p y /2). 

The spectrum of this model is determined by the solutions of the Heisenberg 
equation of motion, ifidt^ pa = [^ P a, Hbk], which is equivalent to solving for the 
LCAO amplitudes: 



(H p - el) * pQ = 0. 



(6) 



This matrix equation is easily solved. The four eigenvalues are determined from the 
secular equation 



det (H p - el) = A(e)xy + B(e) (x + y)+C(e) = 



(7) 
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in which x = s,\T?{p x /2), y = sin 2 (pj,/2) and the coefficients are 

■A(e) = 16 (4:tp d t S p + 2t sp tpp£d — 2tp ( jtpp£ s — tpp£d£s) 
B(e) = -4e p (t%e d + t 2 pd e s ) , 



C(e) 



(8) 
(9) 
(10) 



In the above, the energies £d, £ s and £ p are measured relative to their respective 
atomic levels: £d — e — e^, £ s = e — e s and e p = e — e p . 

To each of the four solutions of equation ([?]) there corresponds an eigenvector 
fp a where B, the band index, labels the solutions. The general (unnormalized) 
form for the eigenvector is 



( D l 


\ 


( 








x i 






\ y } 


J 


\ 



(11) 



4:£pt sptpd y& V ) 

-(£ S £ P - St 2 p y)t pd s x 

8f- tp{l<>y J 

where t 2 p = t 2 p — £ s t pp /2, and e d , £ s and £ p are understood to depend on the 
appropriate value of e, which in turn depends upon p and the model parameters. 
We denote the normalized eigenvector with a tilde: \Pp Q = ^p a /\^p a \- 



2.2. Heitler- London interaction 

The second-quantized Heitlcr-London interaction Hamiltonian describes the in- 
teratomic two-electron exchange between the copper and oxygen orbitalsBEl It 
comprises three parts corresponding to Cu4s <-> 02p a , 02p a <-> Cv&d x a_ y a and 
02p x <-> 02p y exchanges with transition amplitudes J sp , J p d and J pp respectively: 



pp 

HL> 



where, 



■' sp 

~2~ 



E E [^na^n l 3'S , n/3^nc« + ^na^n/J^n^^iic 



(12) 
(13) 



a/3 n 



jrpd 



Jpd 



E E [^L^ljAvfAa + dLy^d^y^ 



(14) 



a/3 n 



^ , (n^ + l,n H )a^n/3^(ri x + l,n„)/3^nQ + ^(na,n v + l)a^n/3^(n«,n s +l)/3^n« 
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E E [^na^n/na + ^L^(! lx + l,,j J^n/3^(ri x + l,n,,)a (I 5 ) 



a/3 n 



I Y"t V't v v 

^ A (n I ,n, + l)a J n^ A (n I ,n s + l)^ J na 



J. 



F. 



+ 1)q (n x + l,n H )/3 (rix,"!,-!- 1 )/ 3 + n„)a 
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The complexity of this expression is greatly simplified by moving to the Bloch basis 
defined in equation (^|). In this momentum-space representation it is more natural 
to label the fermionic operators by the band index, B, than by the orbital index. 
Accordingly we separate these indices by writing 



pa 



/ , pa L pa' 



(16) 



where ^ pa is the normalized LCAO eigenvector coming from the band structure 
and c pa is the annihilation operator for an electron in the -B-band with momentum 
p and spin a. In what follows, we will drop the band-index and focus solely on the 
strongly hybridized <i-band which contains the Fermi level. 

3. BCS Theory 

Moving to the Bloch basis and substituting for ^p,, leads to the momentum-space 
representation of the Heitler-London interaction term. At this point we have made 
no approximations, aside from neglecting completely filled or empty bands, and the 
representation is exact. However, the resulting Hamiltonian is too complicated for 
us to proceed further without approximations. The simplest course we can take 
is to make the mean-field or BCS approximation whereby fluctuations about the 
homogeneous, currentless equilibrium state are ignored. This leads to a simplified 
form of the Hamiltonian: 



(17) 



a/3 pp' 



where the BCS pairing-potential is 



V(p,p') = [J sp S p S p , + J pd D p D p ,\ [W xx ,X p X p , + W yy ,Y p Y p , 

2^ ^-P -^P' W xx ' Wyyl Y p Y p i . 

Here we have introduced a new notation: 



W xx > =c x c' x - s x s' x = 4 cos 



Wyy' = CyCy - S y S' y = 4 COS 



Px+P x 

2 

Py+Py 



(18) 



(19) 
(20) 



Note that Wxx' an d Wyy ar e sums of separable terms, which is to say that the p x 
and p' x dependencies can be factored out. Thus V(p, p') is also a sum of separable 
terms. 

For a sufficiently strong, attractive potential, and a low enough temperature, 
the BCS Hamiltonian gives rise to a superconducting condensate of Cooper pairs 
with zero center-of-mass momentum. Following the BCS prescription we obtain the 
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following self-consistent expression for the superconducting gap:l 



i 



A( P ) 



1 



£*W) 



tanh 




a(p'), 



(21) 



N 



p' 



2E(p>) 



where the quasiparticle energies are -E'(p) = \/£, 2 {p) + |A(p)| 2 , and £(p) = e(p) — /i 
is the normal state energy in the d-band relative to the chemical potential or Fermi 
energy. 

4. Discussion 

In this paper we derived the matrix elements of different interatomic two-electron 
exchange processes and obtained the corresponding gap equation. The analysis 
demonstrates that the B\ g symmetry of the order parameter,! required for the 
description of the ARPES and Josephson data for the gap, cannot be obtained by 
considering only interatomic interactions. We thus arrive at the conclusion that 
other exchange processes should be found and taken into account to explain the 
pairing interaction in the cuprates, if we wish to follow the traditional BCS scheme. 
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